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Abstract 

Results of large-scale Monte Carlo simulations of three-dimensional Ising models with edges and corners are 
reviewed. At the ordinary transition, angle dependent critical exponents are observed, whereas at the surface 
transition edge and corner critical exponents are non-universal and depend on the details of the model. The results 
obtained at the surface transition are compared to exact findings on critical two-dimensional Ising models with 
different types of defect lines. 



1. Introduction 

Critical phenomena at perfect, flat surfaces of 
three-dimensional systems have been the subject 
of intensive research during the last three decades 
[1-3] . For the three-dimensional semi-infinite Ising 
model two different ferromagnetic couplings are 
usually introduced, depending on whether the 
neighbouring spins are both located at the surface, 
Js > 0, or not, Jf, > 0. The resulting phase dia- 
gram is well established. If the ratio of the surface 
coupling Js to the bulk coupling Jf,, r = Js/Jb, 
is sufficiently small, the system undergoes at the 
bulk critical temperature Tc an ordinary transi- 
tion, with the bulk and surface ordering occurring 
at the same temperature. Beyond a critical ratio, 
r > r^p « 1.50 [4,5], the surface orders at a higher 
temperature > Tc at the surface transition, fol- 
lowed by the extraordinary transition of the bulk 
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at Tc. At the critical ratio fsp, one encounters the 
multicritical special transition point. 

The semi-infinite Ising model with a flat surface 
may be considered to be a special case of a more 
complex wedge geometry where two planes meet- 
ing at an angle 9 form an inflnite edge. For 9 — 
77 the flat surface is recovered. Cardy [6] showed 
that at the ordinary transition edge critical expo- 
nents depending continuously on 9 arise on purely 
geometrical grounds. For a given opening angle 9, 
however, the values of the critical exponents are 
expected to be universal and independent on mi- 
croscopic details like the strengths of the coupling 
constants or the lattice type. Whereas edge singu- 
larities at the ordinary transition have been studied 
intensively, especially in two dimensions [7], edge 
critical behaviour at the surface transition of three- 
dimensional Ising models has been largely over- 
looked. At the surface transition the critical fluc- 
tuations are essentially of two-dimensional char- 
acter. One may then argue that the edge should 
act like a local perturbation in a two-dimensional 
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system. In analogy with exact results obtained for 
two-dimensional Ising models with defect lines [8] , 
intriguing nonuniversal edge critical behaviour is 
therefore expected. Similarly, corner critical expo- 
nents depending on the details of the model should 
also be observed at the surface transition. 

In this contribution, I discuss results on edge 
and corner critical behaviour obtained in large- 
scale Monte Carlo simulations of three-dimensional 
Ising models [9-11]. After presenting some details 
of the simulations in the next Section, I then dis- 
cuss edge and corner criticality at the ordinary 
(Section 3) and at the surface transition (Section 
4). A brief summary concludes the paper. 



2. Models with edges cind corners 

To introduce edges in Ising magnets defined on 
simple cubic lattices, periodic bomidary conditions 
along one axis, the 2;-axis, are applied. The remain- 
ing four free surfaces of the crystal may be ori- 
ented in various ways leading to different opening 
angles 9 at the edges. As shown in Figure 1, pairs 
of (100) and (010) surfaces lead to four equivalent 
edges with opening angles 9 = tt/2. The intersec- 
tion of (100) and (110) surfaces form two pairs of 
edges with 9 = ^/4 and 9 = Sir /A, see Figure 2. Be- 
sides the bulk coupling Ji, and the surface coupling 
Js, further couplings may be introduced [9,10]: two 
neighbouring edge spins interact with the edge cou- 
pling Je, whereas an edge spin is coupled to its 
neighbouring surface spin by the edge-surface in- 
teraction Jgs ■ The Hamiltonian is then given by 

bulk surface 

— ^ ^ Jes^rSr' ~ ^ ^ JgSrSr' (1) 

edge— surface edge 

with spins Sr = ±1 at sites r = {xyz). For the 
study of corners we use free boundary conditions in 
all three directions. The couplings between corner 
spins and neighbouring edge spins are set equal to 
Je, thus regarding corners as endpoints of edges. 

The results presented in the following have been 
obtained with the efficient one-cluster-flip algo- 
rithm [12]. Systems with Lx Lx N sites were con- 
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Fig. 1. GcoinctiA' of an Ising model with (100) and (010) 
surfaces forming edges with opening angles d = 7r/2. 
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Fig. 2. Geometry of a model with edges with opening angles 
Q = 7r/4 and 9 = 37r/4. Shown is a cut through the crys- 
tal perpendicular to the edge direction. The full lines rep- 
resent the couplings between neighbouring spins, whereas 
the (110) surfaces are indicated by the dashed lines. 

sidered, with L x L being the number of sites in 
the planes perpendicular to the z-axis, and N be- 
ing the number of sites along that axis. L ranged 
from 10 to 80 and N from 10 to 640. For the in- 
vestigation of corner critical behaviour Ising cubes 
with up to 80^ spins were simulated. 

As we aim at studying properties of systems with 
infinitely long edges, finite-size effects have to be 
monitored closely. The Monte Carlo system should 
be large enough to reproduce the thermodynamic 
values of the bulk and the surface magnetization 
sufficiently far away from the edge. In general, at 
a given temperature, edge quantities must be sta- 
ble against enlargening the system size. The same 
holds for corner quantities. 

An interesting quantity in systems without cor- 
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ners is the magnetization per site for lines parallel 
to the ^;-axis: 

m{x,y) = {\Y,S^yA)/N. (2) 

z 

The sum runs over the N spins in a line, with x 
and y being fixed. The brackets denote thermal 
averages. As visual, the absolute value is taken to 
avoid vanishing magnetizations for finite systems. 

Prom the line magnetizations, various magneti- 
zations of interest are obtained. The edge magne- 
tization TO2 is identical to 

m2 = {\Y.^-ovoz\)IN (3) 

z 

where [xoyoz) denotes an edge site. The surface 
magnetization results from the line magnetization 
in the center of a surface, whereas the bulk mag- 
netization is given by the line magnetization in the 
middle of the crystal. 

The local magnetization mi{x,y,z) may be 
obtained by computing the correlation function 
{Sxyz Sx'y'z') between topologically equivalent 
sites {xyz) and {x'y'z') with maximal separation 
distance. In this way, the corner magnetization of 
Ising cubes with linear dimension L is given by 

ma = V{SoooSlll)- (4) 

I refer to corner quantities by the subscript 3 
(indicating that three planes meet at this point), 
whereas edge quantities have the subscript 2. 



3. Ordinary transition 

In [6], Cardy considered d-dimensional systems 
with {d — l)-dimensional hyperplanes meeting at 
an angle 9. At the bulk critical point local crit- 
ical exponents changing continuously with 9 are 
already obtained in mean-field approximation. 
Cardy showed that all edge critical exponents can 
be obtained by combining bulk and surface crit- 
ical exponents with a new angle dependent edge 
exponent. He also computed the value of this ex- 
ponent in first order of an e = — 4 expansion. 
From the renormalization group point of view, the 
angle dependence has its origin in the invariance 
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Fig. 3. Edge magnetization versus temperature for various 
opening angles 9: tt/4 (circles), 7r/2 (squares), 37r/4 (dia- 
monds), and TT (triangles). Only data not affected by fi- 
nite-size effects are displayed. Error bars are smaller than 
the size of the symbols. 

of the edge under rescaling. This makes the open- 
ing angle to a marginal variable, which may then 
lead to angle dependent local critical exponents. 

Edge magnetizations computed by Monte Carlo 
simulations are plotted in Figure 3 for four different 
opening angles: 9 = 7r/4, 7r/2, 37r/4, and tt (corre- 
ponding to the free surface), with all coupling con- 
stants set equal to Ji,. Only data presenting no no- 
ticeable finite-size dependences are displayed. This 
is achieved by studying systems with a large num- 
ber of edge sites, the edge length being the crucial 
quantity [9]. 

At the ordinary transition, the edge magneti- 
zations vanish at the bulk critical temperature 
ksTc/Jb ~ 4.5115 with a power law behaviour 
m2 ~ {Tc — . Here, is the critical exponent 
of the edge magnetization. In order to estimate 
j32, one may study [13] an effective temperature 
dependent exponent (3'2,eSf{t) defined by 

/32,e//(t) = dlnma/dlnt (5) 

where t =\T — Tc\/Tc is the reduced temperature. 
Certainly, on approach to Tc, t — > 0, P2,eff be- 
comes the asymptotic exponent ^2- Very accurate 
Monte Carlo data are required to get reliable esti- 
mates for the critical exponents. 
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Table 1 

Predictions for the edge magnetization critical exponents 
02{d) from various methods. MF: mean-field approxima- 
tion, RNG: renormalization group theory, HTS; high tem- 
perature series expansions, MC: Monte Carlo simulations. 
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2.50 


1.50 


1.17 


1.00 


RNG [6] 


2.48 


1.39 


1.02 


0.84 


HTS [14] 


2.30 
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0.98 


0.81 


MC 


2.30 it 0.10 


1.28 ±0.04 


0.96 it 0.02 


0.80 it 0.01 



As the slope of /?2,e// changes rather mildly at 
sufficiently small reduced temperatures, meaning- 
ful estimates for the asymptotic exponent /32 {0) are 
feasible. The values for (3'2{0) obtained in this way 
[9] are listed in Table 1 and may be compared with 
results obtained from mean- field theory [6] , renor- 
malizations group calculations to first order in e 
[6] , or high temperature series expansions [14] . The 
predictions of the renormalization group seem to 
be systematically too large, as suggested both by 
high temperature series expansions and our Monte 
Carlo simulations. The last two methods yield re- 
sults which are in close agreement with each other. 

Changing the values of the coupling constants 
or rotating the crystal about the z-axis at fixed 
opening angle 6 (the wedge is then formed by a 
different pair of surfaces) does not alter the criti- 
cal exponents [9] , in accordance with the expected 
universality at the ordinary transition. 

A similar analysis of the corner magnetization ef- 
fective exponent yields for an Ising cube the corner 
magnetization critical exponent /Js = 1.77 ± 0.05 
[11], significantly lower than the mean- field value 
Pi"^ = 2 [7]. 



4. Surface transition 

When approaching the surface transition from 
high temperatures, the surface of a semi-infinite 
system orders whereas the bulk remains disor- 
dered. In a three-dimensional system the critical 
fluctuations are then essentially two-dimensional, 
and the surface critical exponents reflect the re- 
duced dimensionality by taking the values of the 



corresponding two-dimensional bulk system. At 
this transition, an edge may be viewed as a local 
perturbation, acting presumably like a line defect 
in a two-dimensional system. 

Two-dimensional Ising models with defect lines 
have been shown to display non-universal critical 
exponents [8]. A chain defect is formed by a column 
of perturbed couplings with strength Jch, whereas 
at a ladder defect modified couplings of strength 
Ji connect spins belonging to two neighbouring 
columns. In both eases local magnetization critical 
exponent depend on the strengths of these defect 
couplings [8]. In analogy with these results, one 
expects at the surface transition, for a fixed open- 
ing angle, intriguing non- universal edge and corner 
critical behaviour in three-dimensional Ising mod- 
els [10,11]. 

In the following, I discuss the influence of edges 
with opening angle 9 = 7r/2, see Figure 1, on the 
local critical behaviour at the surface transition. 
The ratio Js/Jt ~ 2 > r^p is kept fixed, whereas 
the edge-edge, Je, and the edge-surface, Jes, cou- 
plings are varied. To determine edge critical be- 
haviour accurately, the critical temperature of the 
surface transition, Tg, needs to be determined ac- 
curately. Using standard finite-size analysis, one 
obtains ks Ts/Jb = 4.9575 ± 0.0075 for = 2 Jb 
[10]. Because edges are one-dimensional, and all 
couplings in the models are of short range, edge 
quantities become singular at T,. 

The profile of the line magnetization at the sur- 
face reflects the influence of bulk spins close to the 
surface and edges as well as the strength of the cou- 
plings near the edge, Je and Jea- Typical profiles 
are shown in Figure 4 for Je = Jes = Js = 2J(,. 
The edge magnetization is m(l, J). At low temper- 
atures, m(x, L) increases monotonically with the 
distance from the edge. The magnetization of the 
surface spins, which are connected directly to an 
ordered bulk spin, is enhanced compared to the 
edge magnetization. In contrast, on approach to 
Tg, the ordering of the spins falls off quickly by go- 
ing from the surfaee to the bulk, and the surface 
magnetization is pulled down below the edge mag- 
netization by the coupling to the disordered bulk 
spins. Roughly at Tg, a non-monotonic behaviour 
in the profile shows up, with a maximum close to 
the edge. It may be explained by the fact that sur- 
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Fig. 4. Line magnetization profiles near an edge with open- 
ing angle = n/2 at various temperatures. The couplings 
are Je = Jes = Js = ^Jb- Systems with 40^ spins have 
been simulated. 

face spins next to but on different sides of an edge 
arc more strongly connected to each other, through 
the same neighbouring bulk spin, than spins on the 
flat part of the surface with the separation distance 
of two [10]. 

The effect on the edge critical behaviour of 
weakening the interactions between the edge and 
the surfaces, Jes, is demonstrated in Figure 5. 
A plateau-like behaviour is observed close to Tg, 
which facilitates a precise determination of the 
asymptotic value. Obviously, the value of the crit- 
ical exponent is non-universal, changing continu- 
ously as a function of the coupling strength Jes, 
see Table 2. A similar dependency is found when 
Jes is kept fixed whereas Jg is changed [10]. The 
findings of Figure 5 and Table 2 may be related 
to the reported results on two-dimensional Ising 
models with a ladder defect [8]. The ladder column 
corresponds to the edge, and the ladder couplings 
Ji reflect not only the edge-surface interaction 
Jes but also the reduced connectedness to bulk 
spins at the edge compared to the surfaces. For 
Jes = Je = Js, the critical exponent of the edge 
magnetization has the value (32 = 0.095 ± 0.005, 
significantly lower than the critical exponent of 
the perfect two-dimensionallsing model, /3 = 1/8. 
Of course, the non-monotonic profiles of the line 
magnetization suggest that a closer analogy be- 



Fig. 5. Eflfective edge cxijoiiciits as function of the reduced 
temperature for various values of the coupling Jes with 
Je = Js = "^Jb- Only data not aifected by finite-size effects 
are displayed. The largest system simulated contains 80^ 
spins. Error bars result from averaging over several real- 
izations. 

twecn edge; properties and descriptions by two- 
dimensional Ising models with defect lines would 
require more complicated, extended ladder-type 
defects. 

Corner criticality in three-dimensional Ising 
models deserves to be analyzed at the surface 
transition as well. As edges are local perturbations 
acting similar to defect lines in two-dimensional 
models, the corners of a cube may be interpreted 
as intersection points of three defect lines. 

As for the edge, the local magnetization profiles 
show a non-monotonic behaviour along paths from 
the edges or corners towards the center of the sur- 
face, reflecting the influence of bulk spins. At lower 
temperatures, the profil is monotonic, with the 
smallest magnetization at the corners. Close to Tg, 
a monotonic behaviour is also expected, with the 
highest magnetization at the corners. Due to finite- 
size effects, this last behaviour is not observed in 
the simulated finite cubes [11]. 

The analysis of corner magnetization effective 
exponents yields again critical exponents chang- 
ing continuously with the strengths of the local 
couplings, see Table 2. Note that corner critical- 
ity has not been studied for the case Jes = Js/^- 
Changing corner critical exponents are also ob- 
served when varying Je [11]. This non- universal 
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Tabic 2 

Edge and corner magnetization critical exponents /32 und 
/33 at the surface transition for systems with opening angles 
d = it/2 and Je = Js = 2Ji,. 





Jes — 'Js 


Jes — Js /'^ 


Jes = Js/i 




0.095 ± 0.005 


0.176 ± 0.005 


0.244 ± 0.005 


/33 


0.06 ± 0.01 


0.26 ± 0.02 





corner criticality is understood by the analogy to 
two-dimensional Ising models with star defects 
formed by three intersecting ladder defects [15]. 

5. Conclusion 

Critical phenomena may occur not only in the 
bulk but also at surfaces, edges, and corners. 
I have presented results of an extensive Monte 
Carlo study of edge and corner criticality in three- 
dimensional Ising systems. At the ordinary transi- 
tion, edge and corner critical exponents have been 
computed and compared to analytical predictions. 
Whereas the critical exponents depend on the 
wedge angle, they do not depend on microscopic 
details as the strengths of the local couplings. 
This is different at the surface transition where 
coupling dependent edge and corner critical expo- 
nents are observed. This non-universal behaviour 
is understood by analogy with exact results on 
two-dimensional Ising models with defect lines. 
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